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(54) Determining absorption coefficients or modified scattering coefficients 

(57) A method of determining the absolute absorption coefficient and the modified scattering coefficient in 
highly scattering medium, comprises passing light Into the medium, detecting diffuse light leaving the 
medium, measuring changes in at least two of the intensity, phase and modulation depth that are induced in 
the light by a common cause such as by small absorption coefficient changes or changes in source/detector 
distance. These changes can be Induced by a) changes in the chromophore concentration, b) small changes in 
wavelength, c) changes in blood flow and oxygenation in a limb or in the head of a foetus, or d) changes in the 
distance between light source and detector. The absorption coefficient may be determined from the ratio of 
the two measured changes. 
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De-termlnatlon of absorption coef f iclent:s and modified scattierlng 
coef f Iclentis 

This Invention relates to the determination of absorption 
coefficients and modified scattering coefficients In highly 
scattering media. The Invention finds applications In biomedical 
optics, for example the In vitro or In vivo determination of the 
concentration of various constituents and chromophores (eg. Hb, 
HbOj or oxygen concentrations) In muscle, brain or other tissue. 
It also Is of use In the determination of absorber concentrations 
In other turbid media, for example In food and pigment 
production, or In chemical processes as an aid to manufacturing 
process control. Furthermore, the Invention Is of use for the 
determination of modified scattering coefficients In light 
scattering media, for the determination of concentrations of 
(non-absorbing) constituents like electrolytes, glucose, lactate 
or urea In these media and for the determination of refractive 
Indices . 

The analysis of absorption and scattering properties are an 
Important problem In many areas of biomedical optics. In 
particular for the determination of hemoglobin (Hb), 
oxyhemoglobin (HbO,) and cytochrome oxidase concentrations and 
glucose concentrations In tissue. The different near Infrared 
absorption spectra of these chromophores allow changes In their 
concentrations to be calculated from changes In the intensity of 
light diffusely reflected from the tissue surface. However, the 
absolute concentrations of these chromophores In tissue, le. the 
absolute absorption coefficient, cannot easily be Inferred. 

Different experimental approaches have been suggested for the 
measurement of absorption coefficients (p.) In highly scattering 
media most of which are based upon the measurement of both the 
Intensity and arrival time (time of flight), or phase of an 
Intensity modulated light source, of the light diffusely exiting 
from tissue. It has been shown that the absorption and modified 
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scattering coefficients (p.. P.') can be Inferred from 
roeasurementrs of diffuse Intensity and phase for different source 
detector distances for a fixed modulation frequency. 
Alternatively, single distance, multiple modulation frequency 
measurements can be used. For multiple small source detector 
distances the reflectance measurements are sufficient to derive 
both and p, * . 

In the preferred embodiments of this invention, a different 
approach for the determination of absorption coefficients Is 
adopted. We have found that small changes in absorption 
coefficient or other parameters of a system iinder study (eg the 
distance between a light-source and a detector) induce changes 
in diffuse (reflected or transmitted) light intensity and time 
of flight (or phase «) and modulation depth M of an intensity 
modulated 1 ight wave , and that the rat lo of these changes 1 s 
primarily Independent of the scattering properties of the medium. 
This ratio gives a good estimate of the absolute absorption 
coefficient. This particular approximation is only valid over 
a certain range of and p.* values. However, this range 

encompasses the range found in biological tissues in the near 
Infrared. The term "light" used in this specification therefore 
Includes the infrared (particularly the near Infrared) as well 
as the visible spectrum. The modulation depth M is the ratio of 
the amplitude of the alternating to the steady component of an 
intensity modulated light wave. 

The modified scattering coefficient can readily be inferred from 
the absorption coefficient and from the intensity, time of flight 
and or phase or modulation depth measurement * 

Thus in one aspect the invention provides a method of determining 
an absorption coefficient or a modified scattering coefficient 
of a scattering medium comprising passing light (as herein 
defined) into the medium, detecting diffuse light emitted from 
the medium, measuring changes in at least two of the intensity, 
the phase (or transit time) and the modulation depth of said 
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light: arising from a coxamon cause, and detzermlnlng l:he 
coefficient from said changes* 

In another aspect the Invention provides a method of determining 
an absorption coefficient or a modified scattering coefficient 
of a scattering medltim comprising passing light (as herein 
defined) Into the medium, measuring a change In diffuse light 
emitted from the medium resulting from a small change In the 
absorption coefficient or In the relative positions of a source 
and a detector of the light and determining the coefficient from 
the said change. 

The Invention also provides apparatus adapted and configured for 
use In the methods of the Invention. 

Thus without limitation of the foregoing generality In another 
aspect the Invention provides apparatus for determining an 
absorption coefficient or a modified scattering coefficient of 
a scattering medium comprising means for passing light (as herein 
defined) Into the medium, means for detecting diffuse light 
emitted from the medium, means for causing a change In a 
parameter of the system formed by the apparatus and the medium, 
means for measuring resulting changes in at least two of the 
Intensity, the phase (or transit time) thereof and the modulation 
depth of the emitted light, means for measuring said at least two 
changes, and means for determining said coefficient therefrom. 

In a further aspect the invention provides apparatus for 
determining an absorption coefficient or a modified scattering 
coefficient of a scattering medium comprising means for passing 
light (as herein defined) into the medium, means for inducing a 
small change in the absorption coefficient or for producing a 
small change in the relative positions of a source and a detector 
of the light, means for measuring a resulting change in diffuse 
light emitted from the medium and means for determining the 
required coefficient from said resulting change. 
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Preferably the absorption coefflcieni: Is determined £rom the 
ratio of the change in the emitted diffuse light intensity and 
the change in the transit time (phase) thereof, or from the ratio 
of the change in Intensity or phase to the change in modulation 
depth. 

The small change in absorption coefficient may be induced by 
changing the wavelength of the light passed in to the medium for 
example by a few nanometres. Alternatively the small change in 
absorption coefficient may be induced by varying the 
concentration of a chromophore in the scattering medium. 

The method may also include determining a method as claimed in 
any preceding claim comprising determining the modified 
scattering coefficient by measuring the intensity or phase or 
modulation depth of the input emitted from different locations 
in the medium. 

The determining means may determine the absorption coefficient 
from the ratio of the changes in intensity and transit time. 

There may be means for measuring a phase change in the emitted 
light as the change in traurxsit time. 

There may be means for varying the wavelength of the light by a 
small amount to induce the change in the absorption coefficient. 

In one embodiment the means for passing light into the medium 
comprises a laser diode, the means for varying the wavelength 
comprising means for varying the operating temperature of the 
laser. 

In another embodiment the means for passing light into the medium 
comprises a plurality of light sources having closely adjacent 
wavelengths, the means for varying the wavelength comprising 
means for switching selectively between said light sources. 
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In a furtrher embodiment: t:he means for passing light Intro trhe 
medium comprises a source of whlt:e light:, t:he means for varying 
t:he wavelengt:h comprising means for selec1:lng Individually a 
plurality of closely adjacent wavelengths from the white light. 

The light source and the detector may be a moveable apparatus as 
claimed In claim 11 wherein the light source and/or the detector 
Is moveable and/or there are a plurality of light sources and/or 
detector at spaced apart locations* 

other features of the Invention are as set out In the claims. 

The following Is a more particular but non-llmltlng description 
of the Invention, Including some examples, and refers to the 
accompanying drawings, wherein: 

Figure 1 - shows and as a function of \x^. 

Figure 2 - shows lines of constant as a function of and p.* 
calculated for v„ - 200MHz, n - 1,33 and r = 30mm. 

Figure 3 - shows « ( dA/ar )/0*/ar ) and - 8*/ar )/( dM/ar ) 
as a function of ji. calculated from equations 8 and 9 for r » 
35ram, 0.75, 1.0, 1.25 and l.SOmm'^ and n - 1.56. 

Figure 4 - shows lines of constant as a function of p. and 
* calculated for v„ = 200MHz, n « 1.56 and r = 35mm. 

Figure 5 - a) shows experimentally measured changes In 
attenuation AA, phase A* and modulation depth for variations In 
wavelength. In b) the correlation AA, phase A* Is shown. In c) 
the correlation of A^ and AM Is shown. 

Figure 6 - (a) shows experdLmentally measured changes In 
attenuation AA, phase A*, and modulation depth AM with changes 
In the absorption coefficient p^. In b) the correlation A*, and 
AA Is shown and In c) the correlation of A* and AM. 
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Flgxire 7 - shows absorptrion coefficient = + of the 

liquid phantom, where Is the water absorption and p.** Is the 
dye absorption for a concentration of » 1. 40x10'' v/v. The 
dots are the experimental absorption coefficients derived from 
experiment for the four laser wavelengths. 

Figure 8 - (a) shows the change of attenuation AA and phase 
during a cuff occlusion of the forearm of a volunteer, for optode 
distance r » 40mm, X = 740nm. Figure 8 (b) shows the correlation 
between AA and A* for the data shown in (a) for 50s < t < 170s. 
The slope of a first order regression is AA/AO = 4.67 OD/rad. 

Figure 9 - shows the absorption coefficient p. (solid line) and 
modified scattering coefficient p.' (dashed line) of a solid 
phantom (errors ± 1% and ± 2%, respectively). The absorption 
coefficients derived from Q,(AX) (see Figure 5, filled dot) and 
(see Figure 10, filled squares) are shown. The error bars 
represent the variations by assuming modified scattering 
coefficients between 0.75mm"^< p. '< 2mm"^. 

Figure 10 - (a) shows attenuation changes (AA) and phase changes 
(A*) as a function of source-to-detector distance r measured on 
a solid phantom (p. = 0.0135mm"^, p.'= 0.94mm"^, n = 1.56 for X « 
744nm, compare with Figure 9). The solid lines are the first 
order regression lines through the data and the dashed lines were 
calculated by diffusion theory (equation 8 and 9) and (b) shows 
the coirrelation of AA and A* for the data shown in ( a ) . The first 
order regression has a slope AA/A* =2.38 OD/rad. The absorption 
coefficient p. derived from this slope is shown in Figure 9. 

Figure 11 - shows experimentally determined correlations between 
AA and A*, and between A* and AM, measured on the head of a 
fetus. The data was collected during one contraction. 

Figures 12 and 13 - show various forms of apparatus according to 
the invention. 
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The transport: of light: in 8cat;terlng media may be analyzed by 
diffusion theory to describe light intensity, time of flight and 
phase and modulation depth in terms of the modified scattering 
coefficient the absorption coefficient ji, and the refractive 

index n of the medium. For a "pencil beam" light source on a 
semi- infinite half space, the reflectance R and the mean transit 
time ("time of flight") <t> detected at a distance r from the 
source can be written as 

P 2n-p^ 



and 
[184e 

<t>(p)=- (2) 



2c*(i>+p*^/ii7D) 



respectively. In equation (1) and (2) p=(r^^-2o^)^/^, Zo"l/p,', 
c«Co/n is the velocity of light in the medium (speed of light in 
vacuum Cq), p.„-[3-|i, (P.+P,' )]^^^ is the effective attenuation 

coefficient and D« ~ — — is the diffusion coefficient. 



Here the changes of the reflectance and the mean time with 
respect to changes in are considered whilst the scattering 
properties are assumed to be constant. Any changes in the 
absorption coefficient (Ap.) changes the reflectance from to 
R-Ro+AR(Avi.) and the mean time from <t>o to <t> « <t>o + A<t>. 
The derivative of the attenuation or change in intensity A, 
define by A=log(Ro/R), and <t> with respect to p. are 
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and 

i<^= Zl - (■£♦ -1) (4) 



For the diffusion approximation to apply, p. << p,' and therefore 
P«if " (3|i«P«' )^'* D )• Using these approximations and 

dividing equation (3) by equation (4) gives 
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For t;he optrical properties and source detector distances here 
contemplated, p - Peff/Z >> and therefore 



9j^a In 10-p 



For p > this expression is in good approximation 

independent of p,*, i,e. it is a linear function of p. only. 

When a frequency domain spectrometer ( FDS ) or intensity modulated 
optical spectrometer ( IMOS ) is used to measure chromophore 
concentrations, the phase * of a light wave, intensity modulated 
at the frequency v„, is measured rather than the mean time <t> • 
♦ and <t> are approximately coupled by a simple linear 
relationship at frequencies less than approximately 200 MHz: 

8 



*=-2ifv^'<t> (7) 



The derlva'txve of a't'tenua-tlon and <'t> with respect to the source 
detector distance r Is: 



|7= , / r * ' <2/p+ti.,^) (8) 

ax p2.(i/p+n^^p Inio 



and 



or 2c-(l+PTJi^^p^ 



Alternatively, when a frequency domain spectrometer is used, the 
phase * and the modulation depth M, ie. the ratio of the ac to 
the dc component of the modulated light wave, can be written as: 



*=tx-t:an-M-%-) (10) 



and 



^=ii^*2lS -exp (11) 



where 
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♦r*-*o-sin(e/2) 



(13) 



♦i=to-cos(8/2) 



(14) 



(15) 



and 



(16) 



The equations given above are valid for a homogeneous semi- 
infinite half space. For other geometries, the equations have to 
be modified accordingly. The method described in this preferred 
embodiment of the present invention^ however, is equally able to 
derive p« from these modified equations. 

In Fig. la, ( 3A/ap,)/( 3*/ap,) , the ratio of attenuation 

change over phase change is shown as a function of the 
absorption coefficient p.. The curves in Fig. la were calculated 
from the frequency equation for modified scattering coefficients 
of = 1.0 mm'^ , 1.25 mm"^ and 1 . 50 mm"^ , a source - detector 
distance r=30 mm, a modulation frequency v„ =200MHz and a 
refractive index of n=1.33. It can be seen from Fig. la that a 
good first order approximation is that Q, varies linearly with 
p^. It can also be seen that the influence of p.* on is small. 
Therefore, under the condition of a fixed geometry, constant 
scattering properties and pure absoirption changes, Q, is 
primarily dependent on p^ only. Consequently, measured 
attenuation and phase changes yield a good estimate of the 
absolute (mean) absorption coefficient. The horizontal lines in 
Figure la give the measured values of for four different 
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waveleng1:hs of t:he spectrrometer ±n t:lie experiment: described 
below. 

To furt:her estimate the dependence of on lines of 

constant Q, are shovm In Fig. 2 as a function of p. and y^'. It 
demonstrates that for a change from ji, ' «0.5 mm"^^ to p.* « 
1*5 mm"^, le. about the range of scattering coefficients expected 
for tissue, the corresponding absorption coefficient varies by 
about 10% only. Hence, the exact knowledge of the scattering 
properties Is not required for the prediction of to this 

degree of accuracy* For a given tissue type p,' Is unlikely to 
vary by as much as 0*5 - 1,5 mm"^, leading to a more accurate p. 
estimation* However, the method suggested here assumes that the 
scattering properties of the medium do not change. 

Figure lb shows the ratio V. = 0*/ap.)/( aM/flp.) , the ratio of the 
phase changes over changes In modulation depth, as a function of 
absorption coefficient. The same optical properties of the 
medium were assumed as In figure la. The Influence of the 
different values of p,' on V, Is negligible. To a good 
approximation, V. Is a linear function of p.. Therefore, figure 
lb allows p, to be derived from measured values of V,- 
Equlvalent considerations are valid for the ratio of changes In 
attenuation and changes in modulation depth. 

Equations 1 to 16 assume a matched boundary between medium and 
surrounding. On applying the physically more correct mismatched 
boundary condition for the calculation of reflectance, mean time 
and modulation depth the values for and V, differ only 

slightly from the data in Fig. 1. 

The quotient = OA/dr )/( a*/ar ) calculated from equation 8 and 
9 is shown in figure 3a as a function of p.. A refractive index 
of n -1.56 and a source detector distance of r« 35 mm was 
assumed. Similar to figure 1, the influence of the modified 
scattering coefficient on is small. Consequently, measuring 
the ratio of changes in attenuation and phase with respect to r 
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gives an est±inat:e of p.. The Influence of p,' on the estimated 
values of p^, for constant values of Q^, is shown In figure 4, 
The uncertainties in \x^, due to variations in p,', are of the 
same magnitude as for (compare with figures 1, 2). 

The ratio = { d^/dr) / { dM/dx) is shown in figure 3b. Similar to 
Q^, the influence of p, ' on is small, allowing p, to be derived 
from a measured value of A^/AM. 

Thus^ if em object is probed with light (as defined), the 
concentrations of the constituent substances within the body can 
be estimated from measurements of the attenuation and pathlength 
(time of flight) or from measurements of changes in pathlength 
and modulation depth or changes in attenuation and modulation 
depth. There are two technologies that may be used to measure 
both the attenuation and pathlength of light through an object* 
As the two methodologies, time domain or frequency domain, are 
distinct each method will be described separately. 

Time Domain Measurement 

The attenuation and pathlength of light through an object can be 
measured most accurately in the time domain by using a picosecond 
laser and a light detection device such as a streak camera. The 
laser is used to deliver an ultrashort light pulse to the object 
and the light detection device is then used to measure how the 
shape of this pulse is changed by its interaction with the object 
(i.e. the detection device measures the impulse response ) . 

If the measured impulse response is a function g(t), with a mean 
<t> given by. 
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f g{t)t dt 



<t>=— 



(17) 



fgit) dt 



the pa-bhlength (or mean time of flight) can then be calculated 
from the expression 



where c^ is the velocity of light in vacuum and n is the mean 
refractive index of the object being probed. 

Systems of the above type are large and expensive and this is 
likely to preclude their use in all but the laboratory setting, 
or as fixed equipment in manufacturing or processing plant. 

Suitable light sources for use in such systems are ion laser 
pximped dye lasers and ion laser pumped solid state lasers. 
Suitable detectors are streak cameras, microchannel plate 
photomultiplier tubes and time resolved single photon counting 
systems . 

Alternatively, because only the mean pathlength or transit time 
need be investigated and not the complete Impulse response of the 
object being probed, it is possible to use simpler 
instrumentation, with a poorer time resolution* This alternative 
method consists of probing the object with a broad pulse • 



(18) 



n 
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Frequency Domain Measuremen't 



To measure the attenuation and pathlength of light through an 
object, via the frequency domain, a modified Intensity modulated 
optical spectrometer (IMOS) is used. With an IMOS the intensity 
of the light is sinusoidally modulated at high frequencies (i.e, 
typically > 1 MHz ) . A phase sensitive detector then determines 
the difference between the phase of the optical waveform entering 
the object and the phase of the optical waveform leaving the 
object. As the measured phase difference is a linear function 
of the pathlength, the pathlength can be easily determined. 

Thus, using either an instrximent for time domain measurements or 
for frequency domain measurements to detect intensity and time 
of flight (or phase) and modulation depth changes, it is possible 
to determine the absorption coefficient and the modified 
scattering coefficient. 

One approach to obtain Q, is to induce small changes in 
absorption coefficient by tuning the wavelength X ie. to scan 
over the absorption spectrum of the scattering medium. It is 
necessary that the scattering coefficient does not change 
materially for these wavelengths. 

To demonstrate this technique, a solid, light scattering phantom 
with well characterised optical properties (for X « 758 nm; p. = 
0*0160 mm; p.' « 0.934 mm; n=1.56; 3^ = +1.5% per nm; ap,'» 
-0.05% per nm) was used. A temperature tuneable diode laser was 
employed to change the wavelength between X = 753 nm and X = 761 
nm (FWHM « 2 nm). The light reflected from the phantom was 
detected at a distance r * 30 mm from the light delivery fibre 
and attenuation and phase changes were corrected for 
instabilities in the diode laser output by reference measurements 
at a distance r = 7mra, 

Figure 5a shows the measured changes in attenuation, phase and 
nicdulation depth for this wavelength range. The first order 
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regression of 1:he correlation plots gives slopes of = 2,23 
OD/rad (figure 5b) and = 5.03 rad (figure 5c). By analysing 
these values according to equations given above, an absorption 
coefficient of p. « 0.0161 min"^ for 0.75 min"^ < p,'<2.0 mm"^ was 
calculated from and = 0.0163 mm'^ from V.. This Is In 
excellent agreement with the true value of p.. 

The wavelengths have to be adjusted to the absorption spectrum 
of the chromophores. The absorption coefficients of hemoglobin 
and oxyhemoglobin for example change by up to * ±2% per 

AX»lmm, I.e. a variation of the wavelength by a few nm Is 
sufficient to Induce large enough absorption changes. 

Thus by modifying the IMOS to permit tuning of the wavelength of 
the light source over a narrow range eg. a few nanometres, for 
example (between 5 and 20 nm) the resulting change In attenuation 
AA, phase A« and modulation depth AM can be measured, enabling 
the absolute absorption coefficient to be determined. Because 
Is determined from AA, A* and modulation depth AM for the same 
AX, the exact value of AX, Is not critical. 

The IMOS light source may for example be a modified laser diode, 
a modulated light emitting diode or a lamp whose output Is 
controllable with the necessary speed and precision. Modulation 
may be by any single or multlf reguency slnusoldally based scheme. 
Suitable detectors for the IMOS are photomultlpller tubes, 
photodlodes and CCD (charge coupled device) detectors. 

The wavelength tuning capability may be achieved in a number of 
ways, for example as already stated by adjusting the temperature 
of the light source, or by adjusting the current supply, or by 
a tunable absorbance filter. Alternatively, a white light source 
may be used with a Pockel cell and other wavelength selective 
device (eg. a prism or grating) as a number of discrete 
selectable light sources with closely spaced wavelengths (e.g. 
multiple LEDs) may be employed. 
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In an experiment: to det:ennlne an absorption coefficient by 
concen-tratlon adjustment, a frequency domain spectrometer was 
used Incorporating four different laser diodes (X«744 nm, 807 nm, 
832nm, 860 nm) that were Intensity modulated with a frequency of 
v„=200 MHz. Phase sensitive amplifiers detected the phase shift 
of the multiple scattered light. Optical fibres were used to 
transport the light between light source, and the scattering 
medium and between the scattering medium and detector. 

Scattering phantoms of known absorption and modified scattering 
coefficient (p, and p,', respectively) were employed to 
demonstrate the method. The phantoms consisted of spherical 
polystyrene particles suspended In water serving as light 
scattering centres. Mle theory was used to derive the modified 
scattering coefficient p.* for these spheres (diameters between 
0.6 ym and 2.5 pm). The absorption coefficient of the phantom 
(p.) was that of water (p/) and of a dye (p/) (S109564, ICI, 
Manchester, UK) added In known quantities: p. = p^" + p.**. The 
spectrtim of p^** + p/ for a dye concentration of c^ = 1.40x10"' v/v 
Is shovm In Fig. 7. The ends of the light delivering and 
detecting fibres of the frequency domain spectrometer were 
submerged approximately 2 mm Into the phantom of a volume V = 100 
mm X 80 mm x 60 mm. 

The absoi*ptlon coefficient p^ of the phantom was changed by 
varying the concentration of the dye In 14 steps of 1.44 % 
each (corresponding to Ap. = 1.56 xlO'* mm"^ at X « 744nm). The 
concentration of the scattering microspheres c, was 1 % v/v and 
held constant. The resulting optical properties of the phantom 
are listed in table 1 for the different laser wavelengths. The 
wavelength dependence of p,' is Insignificant while p. varies by 
about 30 %. The values of p, given in table 1 are for a dye 
concentration of c^ = 1.40»10"* v/v. 
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Optica properties of the phantom for four different 
wavelengths X. The value of p.* was derived from Mie 
theory. p, consists of the dye absorption and the water 
absorption p/. The last two columns give the measured 
ratio of the attenuation and phase change (Aa/A*) and the 
corresponding absorption coefficient derived from Fig. 1. 

The changes in intensi-ty phase and modulation depth as a function 
of changes in the absorption coefficient were recorded using the 
frequency domain spectrometer. This is shovm in Fig. 6a for X 
- 744 nm. On plotting AA against A* (see Fig. 6b) it can be seen 
that both are strongly correlated. A first order regression 
gives a slope of AA/A* - 2.47 OD/rad. By using this measured 
value of Q, = AA/A* and the known modified scattering 
coefficient. Fig. la allows the absorption coefficient y. to be 
read. The experimental values of are listed in table 1 and 
shown in Fig. la for all four wavelengths. The corresponding 
absorption coefficients were read from Fig. la and are plotted 
in the absorption spectrum of Fig. 7. They agree with the 
expected values to within a few percent (compare table 1). The 
correlation between AM and A* is shown in figure 6c. A first 
order regression gives a slope of A*/Am = 3.53 rad. Using this 
value allows to be read from figure lb to be 0.0135 mm"^. This 
is in agreement with the true ji, value as well. 

In a second experiment the method was used to determine an 
absorption coefficient in vivo and thereby demonstrate the change 
of blood flow during venous occlusion. 

Fig. 8a shows attenuation and phase changes for a wavelength 
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X « 740 mm during a cuff occlusion of the forearm of a healthy 
volunteer. The optode spacing was r=40 mm and the modulation 
frequency of the frequency domain spectrometer v„«200 MHz. The 
pressure In the cuff was Increased to about 200 mm Hg at t»60 s 
and released at t=I80 s. Both the attenuation and the phase 
difference increase during the occlusion. Indicating an Increase 
In the absorption coefficient due to a change In the oxygenation 
of the blood. The correlation between attenuation and phase 
change for times between 50 s and 170 s (see Fig. 8b} gives a 
slope AA/A* = 4.85 OD/rad. Assuming a refractive Index of the 
tissue of n « 1.4 equations 3, 4 and 7 give values for the 
absorption coefficient of = 0.026 mm"^ (for * « 0.5 mm'^) and 
« 0.028 mm*^ (for * = 1.0 ram"^)* 

In a further experiment the attenuation and phase changes were 
measured on the head of a fetus during labour. Contractions 
change the blood flow and volume In the baby's head. I.e. Induce 
changes In absorption coefficient. Figures 11a and b show the 
correlation of attenuation and phase changes and the correlation 
of phase and modulation depth changes measured during one 
contraction of about 1 . 5 minutes duration for two wavelengths ( \ 
«= 807 nra and X, « 832 nm). The measured slope AA/A* (Xt.) * 4.72 
OD/rad corresponds to p. (X^^) «= 0.026 mm"^ and the slope AA/A* (Xj) 
= 6.1 OD/rad corresponds to (X^) = 0.033 mm*^. Using the 
absorption coefficients of pure hemoglobin and oxyhemoglobin at 
both wavelengths, the oxygen saturation of the blood of the 
tissue probed by the light can be calculated to be 77%. The 
measured slope of phase and modulation depth changes A*/ AM (X^) 
« 2.75 rad (see figure lib) corresponds to (X^^) = 0.012 mm"^. 
This absorption value Is different from that derived from AA/A*. 
That Is due to the fact that attenuation and modulation depth 
probe different volumes of the tissue. While attenuation 
measurements are more susceptible to changes In the volume close 
to the light source detector, le. to the surface of the medium, 
the phase and modulation depth probe predominantly volvime deep 
within the medium. A*/ AM Is therefore more likely to give the 
absorption coefficient (and oxygen saturation) of the deeper 
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layers of t:lie brain. Contrary, AA/A* derives the absorption 
coefficient of the outer tissue layers. 

Experiments on phantoms were undertaken to derive p. from 
measurements with varying source detector separations. The 
attenuation and phase changes were measured as a function of the 
source detector distance for the solid scattering phantom (see 
Figure 10a). The optical properties of the phemtom were n - 
1.56, P« - 0.0134 mm"^ and p,' « 0.094 mm"^ for X « 744 nm. In 
figure 10b, the correlation of AA and A4» » 2.38 OD/rad. Analysis 
of this slope according to equation 8 and 9 (compare figure 3 and 
4) yields an absorption coefficient 0.0135 (± 0.0010) mm"^ for 
assumed modified scattering coefficients between 0.75 mm*^ < p.* 
< 2 mm~^. This is in excellent agreement with the true value of 
0.0134 mm*^. 

Using the estimated value of p^, together with the measured AA, 
allows the modified scattering coefficient to be calculated from 
equation 1. The slope of the first order regression, AA/Ar, in 
figure 10a is 0.109 CD/mm (X = 744 nm) giving 0.900 (± 

0.0065) mm'^ for = 0.0135 (±0.001) mm"^. These values are in 
agreement, to within the error limits, with the true p,' . 
Equivalently, the slope of the phase A*/Ar - -0.0460 rad/mm. 
Using this value together with the estimated p^ and eq[uation 2 a 
modified scattering coefficient of p,'« 0.902 (±0.065) mm*^ is 
obtained. 

Figures 12A, 12B, and 12C show diagrammatically various forms of 
apparatus of the invention. The apparatuses of Figures 12B and 
12C are the same as that of Figure 12A except as illustrated. 

In Figure 12A, a sample 10 is supplied with light as hereinbefore 
described via an optical fibre 12 from a source 14. The source 
is driven by control circuitry 16 of a spectrometer 18 so as to 
have known amplitude and phase; the circuitry 16 controls the 
operating temperature or driving current to adjust the wavelength 
of the light over a small range. A detector 20 receives diffused 



19 



light from the sample 10 via optical fibre 22, enabling the 
spectrometer and computing equipment 24 connected thereto to 
determine p. as already described. 

Shown diagrammatically at 25 are means such as a controllable 
supply of solvent for adjusting the concentration of the sample 
as an alternative means of varying the absorption coefficient p^. 

In Figure 12B the light source 14 is replaced by four separate 
sources 26, connected via a coupler 28 to the optical fibre 12. 
The control circuitry 16 switches selectively between them to 
vary the wavelength of the light supplied to the sample. 

In Figure 12C, the light source 14 is replaced by a white light 
source 30, in connection with a Pockel cell, prism or grating 32 
controlled from the control circuitry 16 to provide monochromatic 
light of known and controllable wavelength. 

Figure 13 shows variations in the source/detector configuration 
when a change in the relative positions thereof is employed to 
generate the changes in A, * or M. Except as now described the 
apparatus is as in figure 12A except that the wavelength of the 
light source is fixed. 

In figure 13A, the sample 10 is illuminated by one of four 
selectively switchable light sources 34, the diffuse light being 
received in transmittance (through the sample) by detector 36. 
Selecting different ones of the light source changes by a small 
accoxxnt the distance between the light source and the detector. 
It is not necessary to measure the exact change in distance; it 
is sufficient that it is the common cause of a change in A, * and 
M. 

The apparatus of figure 13B is similar to that of figure 13A 
except that there is a single light source 34 and a plurality of 
sui t ab le detectors 36. S imi 1 ar changes in source - detector 

distance are achieved. 
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In £lg\ire 13C -there Is a single light: source 34 and a single 
detector 36 which is arranged t:o receive diffuse light: re£lect:ed 
out of the same surface of the sample 10 as it entered* The 
source-detector distance is adjusted by moving one or both of the 
source and detector* Alternatively several switchable sources 
or several switchable detectors can be provided as in figures ISA 
and 13 B. 

The main features and advantages of the specifically-described 
embodiments of the invention are: 

A single source detector distance is sufficient to determine 
when the changes in AA , and AM are induced by absorption 

changes are used. The absolute change in distance does not need 
to be known to determine ji, when the changes in AA, A* and AM 
are induced by variations in the source detector distance. 
A single modulation frequency is sufficient. 

Changes in intensity, phase and modulation depth are measured 
rather than absolute values. 

Therefore, the technology needed is simpler than that for other 
methods . 

Each feature disclosed in this specification (which term 
includes the claims ) and/or shown in the drawings may be 
incorporated in the invention independently of other disclosed 
and/or illustrated features. 

The following abstract was filed with the application and is 
repeated here as part of the specification: a method of 
determining the absolute absorption coefficient (chromophore 
concentration) and the modified scattering coefficient in highly 
scattering media, comprising measuring the intensity and phase 
and modulation depth changes that are induced by small absorption 
coefficient changes or changes in source detector distance. 
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CLAIMS 

1. A method of determining an absorption coefficient or a 
modified scattering coefficient of a scattering medium comprising 
passing light (as herein defined) into the meditom, detecting 
diffuse light emitted from the medium, measuring changes in at 
least two of the intensity, the phase (or transit time) and the 
modulation depth of said light arising from a common cause, and 
determining the coefficient from said changes. 

2* A method as claimed in claim 1 wherein the common cause is 
a small natural or Induced change in the absorption coefficient 
or a small change in the relative positions of a light source 
to a detector thereof, 

3. A method of determining an c^Dsorption coefficient or a 
modified scattering coefficient of a scattering medium comprising 
passing light (as herein defined) into the medium, measuring a 
change in diffuse light: emlt:ted from the medium resulting from 
a small change in the absoirptlon coefficient or in the relative 
positions of a source and a detector of the light, and 
determining the coefficient from the said change. 

4. A method as claimed in any preceding claim wherein the 
absorption coefficient is determined from a ratio of two said 
changes . 

5. A method as claimed in claim 1 where the change in the 
transit time is measured by determining a change in phase of the 
emitted light. 

6. A method as claimed in any of claims 2 to 5 wherein the small 
change in absorption coefficient is induced by changing the 
wavelength of the light passed into the medium. 

7. A method as claimed in claim 6 wherein the wavelength is 
changed by a few nanometres. 
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8. A metzhod as claimed In any of claims 2 -to 5 wherein the 
small change In absorption coefficient Is Induced by varying the 
concentration of chromophores In the scattering medium. 

9« A method as claimed In claim 2 wherein the change In 
pathlength or phase Is achieved by moving the light source and 
or the detector. 

10. A method as claimed In any preceding claim comprising 
determining the modified scattering coefficient by measuring the 
Intensity or phase or modulation depth of the light emitted from 
different locations In the medium. 

11 • Apparatus for detezmlnlng an absorption coefficient or a 
scattering coefficient of a scattering medium, said apparatus 
being adapted and configured for use In the method of any 
preceding claims. 

12. Apparatus for determining an absorption coefficient or a 
modified scattering coefficient of scattering medium comprising 
means for passing light (as herein defined) Into the medium, 
means for detecting diffuse light emitted from the medium, means 
for causing a change In a parameter of the system formed by the 
apparatus and the medium, means for measuring resulting changes 
In at least two of the Intensity, the phase (or transit time) 
thereof and the modulation depth of the emitted light, means for 
measuring said at least two changes, and means for determining 
said coefficient therefrom. 

13. Apparatus for determining an absorption coefficient or a 
modified scattering coefficient of a scattering medium comprising 
means for passing light (as herein defined) into the medium, 
means for Inducing a small change In the absorption coefficient 
or for producing a small change in the relative positions of a 
source and a detector of the light, means for measuring a 
resulting change in diffuse light emitted from the medium and 
means for determining the required coefficient from said 
resulting change. 
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14 • Apparatus as claimed in claim 12 or claim 13 wherein -the 
determining means determines the required coefficient from the 
ratio of two said resulting changes. 

15. Apparatus as claim in claim 14 wherein the determining means 
determines the absorption coefficient from the ratio of the 
changes in intensity and transit time, 

16. Apparatus as claim in claim 4 comprising means for varying 
tiie wavelength of the light by a small amount to induce the 
change in the absorption coefficient. 

17. Apparatus as claimed in claim 16 wherein the means for 
passing light into the medium comprises a laser diode, the means 
for varying the wavelength comprising means for varying the 
operating temperature of the laser. 

IB . Apparatus as claimed in claim 16 wherein the means for 
passing light into the medium comprises a plurality of light 
sources having closely ad j acent wavelengths , the means for 
varying the wavelength comprising means for switching selectively 
between said light sources. 

19. Apparatus as claimed in claim 16 wherein the means for 
passing light into the medium comprises a source of white lights 
the means for varying the wavelength comprising means for 
selecting individually a plurality of closely adjacent 
wavelengths from the white light. 

20. Apparatus as claimed in claim 11 wherein the light source 
and/or the detector is moveable and/or there are a plurality of 
light sources and/or detector at spaced apart locations. 

21. A method or an apparatus for determining an absorption or 
coefficient or a modified scattering coefficient substantially 
as herein described with reference to the accompanying drawings - 
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